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Abstract. We investigate the recent claim by Qin et al. that the 
observed correlations (or lack thereof) between the core dom¬ 
inance parameter and the core and extended powers of sam¬ 
ples of lobe- and core-dominated quasars is in contradiction 
with beaming models. Contrary to their conclusion, we find 
that their results are in perfect agreement with such models, 
and support this assertion with Monte Carlo simulations. 

Key words: galaxies; active - quasars: general - radio contin¬ 
uum: galaxies 


1. Introduction 


It is widely believed that radio-loud quasars and powerful radio 
gal axies ( FR Il-type according to the scheme of Fanaroff & Ri¬ 
ley 1974 1 differ only by the orientation of their radio axes w ith 
respect to the observer’s line of sight (e.g. Barthel 1989 ). A 
related scenario proposes that the cores of these radio sources 
are dominated by beamed emission which can be significantly 
enhanced by Doppler boosting, while the extended emission is 
isotropic (e.g. Orr & Browne 1982). These two models explain 
why radio loud quasars tend to have more luminous radio cores 
compared to FR II radio galaxies with similar extended radio 
luminosities, and has prompted the use of the core dominance 
parameter, R = P§/P^ as a relative measure of orientation. 
[In an attempt to have a meaningful, concise, and consistent 
notation throughout this paper, we employ the superscripts c, 
e, and t to refer to the core, extended, and total radio emission 
(power, P, or flux, S) respectively, and subscripts to indicate 
the rest-frame frequency (in GHz).] 


Recently, Qin et al. ( 1998 hereafter Q98) investigated the 
core and extended radio properties of four samples of quasars 
from the literature. They found that in core-dominated quasars, 
the extended radio power, P®, was inversely correlated with R, 
while the core radio power, P®, was not correlated at all. In 
lobe-dominated quasars, on the other hand, P® was uncorre¬ 
lated while P® was correlated with R. They claimed that these 
correlations were in contradiction with the unified schemes and 
beaming models, and proposed instead that the core emission 
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in core-dominated quasars is unbeamed and isotropic. In this 
Letter, we show that the results found by Q98 are in agreement 
with beaming models, given the selection criteria of the sam¬ 
ples they studied, and that core-dominated quasars do not com¬ 
prise a separate class of active galactic nucleus. In §|| we ex¬ 
plain the reasons for the correlations observed (or unobserved) 
by Q98, and in §|| we compare the predictions of Monte Carlo 
simulations of extragalactic radio sources with the data from 
the samples analyzed by Q98. 


2. Beaming models 

In beaming models, it is assumed that all radio sources 
are intrinsically steep-spectrum, and flat-spectrum sources are 
merely those where the core is boosted sufficiently strongly to 
dominate the total flux over the observed frequency range. The 
extended radio luminosity of all radio sources therefore fol¬ 
lows the radio luminosity function for steep-spectrum sources, 
in which the number of sources above a given luminosity de¬ 
creases sharply with luminosity. 


2.7. Core-dominated quasars 

We first consider core-dominated quasars (CDQs). Q98 believe 
that P® should be correlated with R because those sources with 
the most luminous cores should be those with the strongest 
beaming. This is obviously true at some level, since both Pun- 
sly’s (1995) and Murphy, Browne & Perley’s (1993) samples, 
which are selected on core strength, contain only strongly- 
beamed sources. However, we show by considering the situ¬ 
ation in more detail that the presence or absence of correlations 
is controlled by selection effects, and that this simplistic belief 
is too naive. 

Beaming models predict that the observed core dominance 
parameter should depend on the jet speed, /3c, the viewing an¬ 
gle, 9, and the transverse (unbeamed) core dominance parame¬ 
ter, Rt, as 


P = lprp(l —/3cos6>) ^-b 1 pt( 1 +/3cos0) ^ 

(e.g. Orr & Browne |1982 ), where the second term can be ig¬ 
nored for small values of 6. Obviously, the orientation of any 
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given source is independent of its extended radio luminosity, 
and let us assume for the moment that f3 and i?T are also in¬ 
dependent of P®. In this case, the extended radio luminosity 
function (ERLF) is separable in some form p = f{P, z)g{R) 
and sources with a given value of R follow the same ERLF as 
the general source population. Let us examine Punsly’s ( 1995| ) 
sample, which is Q98’s Sample 1. This consists of RLQs with 
core powers (integrated over 10 MHz < v < 250 GHz) above 
1039 W. Consider those sources with, say R 10^; they need 
to have P® > 10^® W to satisfy the selection criterion. The 
shape of the ERLF is such that about three-quarters of all such 
sources will lie within a factor two of this limit, and only about 
1% will exceed it by an order of magnitude. Most sources with 
ii « 10^ will therefore have P® w 10^® W. If one now con¬ 
siders sources with R 10^, the same arguments indicate that 
these sources will have P® « 10^^ W. Therefore the extended 
radio power will be strongly inversely correlated with R, which 
is exactly what is observed in Punsly’s sample. Since the core 
luminosity is merely P® = PP®, the core power will be uncor¬ 
related with R. Exactly the same argument applies for Sample 2 
(from Murphy et al. 1992), which is a sample selected on core 
flux, rather than luminosity. It is the fact that these samples are 
selected on core strength that causes the observed correlations, 
and not the fact that the sources are CDQs. 

The maximum observable core dominance parameter is 
Rmax ~ ^Pt (1 For typical values of/3 = 0.99 (T = 7, 


(S§ > 0.1 Jy), which means that it will include many CDQs. 
In addition, the sample is drawn from the Hewitt & Burbidge 


(1987) optical quasar catalogue and is therefore subject to un¬ 


certain selection effects. The fact that in reality it is a mix of 
CDQs and LDQs can serve to explain why it displays both cor¬ 
relations, but at a much weaker level than in samples of CDQs 
or LDQs alone. We do not consider it further. 

Instead, we turn to the well-defined Sample 4, from Hough 
& Readhead ( 1989| ) which contains all double-lobed quasars 
from 3CRR (Laing, Riley & Longair 19831 whose cores were 


bright enough to be mapped with the Mark III VLBI system. 
All sources therefore have Sq > 10.9 Jy and S§ > 5mJy. 
Unlike the core-selected samples, the total flux at the selec¬ 
tion frequency (178 MHz) has a negligible contribution from 
the flat-spectrum core, and therefore Sq j^^g has no dependence 
on P. Neither does Pg j^^g, which can be assumed to be the 
extended luminosity alone. The spectral index of the lobes is 
not expected to be physically related to core dominance, and 
we confirm that t his is s o in the sample of radio galaxies from 
Hardcastle et al. ( 1998 ) by using the generalized Spearman’s 
rank order correlation coefficient (Isobe, Feigelson & Nelson 


1986), where we And that the two quantities are uncorrelated 


e.g. Guerra & Daly 1997 ) and log Pt = —2.5 (at 5 GHz, Simp¬ 
son 1996; see also Morganti et al. 1997), this is about 20 using 


the defintion and canonical radio spectra of Punsly ( 1995 ). To 
exceed this value by nearly two orders of magnitude (as some 
of Punsly’s sources do), a source must either have an unusually 
fast jet (/3 = 0.999, T = 22) or a much larger than average 
Pt. These parameters will therefore be correlated with the ob¬ 
served core dominance, and so could an additional correlation 
with extended radio power nullify the above argument? 

Suppose that /3 is correlated with the extended radio lumi¬ 
nosity. This is not unreasonable since the most powerful radio 
sources might also produce the fastest jets. If so, sources with 
large P are more likely to be drawn from the high-luminosity 
end of the ERLF. On the other hand, Pt might be anticorre¬ 
lated with extended radio luminosity, since the lobe luminosity 
of a radio source decreases thro ughou t its lifetime despite a 
constant jet power (e.g. Baldwin 1982 ). So sources with large 
P might preferentially be older sources with lower extended 
radio luminosities. The fact that the anticorrelation between 
log P® and log P appears not to alter its slope at large P in¬ 
dicates that either both these effects are weak, or that they have 
approximately the same strength and cancel each other out. 

2.2. Lobe-dominated quasars 

We now turn to lobe-dominated quasars (LDQs). Hooimeyer et 


at greater than 90% confidence. Therefore, P§ will be uncorre¬ 
lated with P. On the other hand, Pq = R P|, and should there¬ 
fore be clearly correlated with P, with the linear regression line 
having a slope of unity, in line with observations. Again, it is 
not the nature of the sources (LDQs) which produces the result, 
but the selection criteria. 


3. Monte Carlo simulations 

In order to test whether the above arguments can explain the 
observations in a quantitative manner, we perform simulations 
of radio sources. We start with the steep spectru m rad io lumi¬ 
nosity function (RLFl) of Dunlop & Peacock ( 1990| ), which 
we assume to be equivalent to the extended luminosity of all 
radio sources, and consider only those sources with luminosi¬ 
ties abo ve the Fanaroff-Riley break (i.e. FRII sources; Fanaroff 
& Riley 1974 ), since it is only these which con stitute the par¬ 
ent population of quasars (e.g. Urry & Padovani 1995 ). We as¬ 
sume that the sources are randomly oriented with respect to the 
observer, and that their log Pt values are drawn from a Nor¬ 
mal distribution with mean —2.43 and standard deviation 0.50. 
This distribution provides a good fit to the observed values of 
P for FRII radio galaxies and quasars from Laing et al. (1983) 
with z < 0.43, after applying the Orr & Browne (1982) beam¬ 
ing model (Simpson et al. 1998; see also Simpson 1996). It 


al.’s (1992) sample of quasars (Q98’s Sample 3) is selected on a 


minimum angular size, thereby requiring measurable emission 
from the radio lobes and excluding the most strongly-boosted, 
pole-on sources, but also has a bright core radio flux limit 


is assumed that the radio cores have a = 0 and the extended 
emission has a = 0.7 (Si, oc A jet speed [3 — 0.99 is 

assumed for the purposes of simulating Doppler boosting. We 
then subject the ensemble of sources to the selection criteria of 
the samples studied by Q98, and compare our regression and 
correlation analysis with theirs. As the samples of Q98 invari¬ 
ably have a number of sources with incomplete information 
that prevent them being included in Q98’s analysis, we also 
randomly exclude sources from our complete samples so as to 
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have the same sample size. This ensures that the same level 
of randomness should exist in our simulated samples as in the 
real data. All our simulations assume Hq = 50kms“^ Mpc“^, 
go = 0.5 and A = 0, although the significance of the correla¬ 
tions we find do not depend on our choice of cosmology. 


Table 1. Correlation coefficients, r, and slopes of the best-fit 
regression lines, a, for the three samples modelled. For each 
sample, the top line contains the results from the real data de¬ 
termined by Q98, and the bottom line contains the results from 
the simulated data. 


3.1. Sample 1 


Sample 1 (Punsly 1995) consists of quasars with core lumi¬ 
nosities, integrated over the range 10 MHz < v < 250 GHz, 
greater than 10^® W. Adopting the typical core spectrum as¬ 
sumed by Punsly, this corresponds to > 8 X 10^’^WHz 
which we adopt as our selection criterion. In Figure we plot 
the extended and core luminosities against the core dominance 
parameter. 


3.2. Sample 2 


Sample 2 (Murphy et al. 1992) contains sources with core 
fluxes iSg > 1 Jy and spectral indices (between 1.4 and 5 GHz) 
a < 0.5 which are identified with quasars. This last criterion 
will exclude nearby sources since they would be classified as 
(broad line) radio galaxies. Examining Murphy et al.’s sample, 
we elect to discard all sources from our dataset with z < 0.25 
to simulate this. There are in fact only about three sources at 
such low redshifts in a typical dataset, so the results are not 
strongly affected. The simulated datasets are shown in Figure^. 


3.3. Sample 3 


As discussed earlier. Sample 3 (Hooimeyeret al. 1992) has un¬ 
certain selection criteria, and so we make no attempt to simu¬ 
late it. 


3.4. Sample 4 

Sample 4 (Hough & Readhead 1989| ) consists of sources with 
Pg > 10.9 Jy and iSg > 5 mJy that are classified as quasars. 
Since this is a much fainter core flux limit than Sample 2, we 
cannot just apply a redshift cutoff, since many distant narrow 
line radio galaxies will have sufficiently bright cores. Instead 
we require (in addition to z > 0.25), that 0 < 45° (see Barthel 
1989), which produces the sample displayed in Figure]^. 


Sample 

N 


rpc^R 

ape^Ji 

apc.R 

1 

94 

-0.899 

0.196 

-0.91 ±0.05 

0.09 ± 0.05 



-0.864 

0.257 

-0.87 ±0.10 

0.13 ±0.10 

2 

54 

-0.739 

0.150 

-0.88 ±0.11 

0.12 ±0.11 



-0.591 

0.136 

-0.84 ±0.10 

0.16 ±0.10 

4 

28 

-0.286 

0.848 

-0.16 ±0.10 

0.84 ±0.10 



0.056 

0.873 

0.03 ±0.10 

1.03 ±0.10 


corporated a range of jet speeds into our model, since with a 
single speed {fS = 0.99) such sources require a viewing angle 
close to the line of sight and a large value of R^. If the jet speed 
is a third random variable then such sources can be created if 
only two out of the three variables have extreme values, and 
this would therefore increase their frequency in our simulated 
Universe. 


4. Summary 


We have investigated the claim by Qin et al. (1998) that 
the observed correlations between core dominance and core 
or extended radio luminosity in samples of core- and lobe- 
dominated quasars are in contradiction with the beaming mod¬ 
els which form a part of the unified scheme for extragalactic 
radio sources. We have argued that these correlations are ac¬ 
tually expected to arise from the beaming models as a natural 
consequence of the steepness of the (extended) radio luminos¬ 
ity function and the selection criteria of the samples. We have 
supported our claim by simulating quasar samples using the ob¬ 
served radio luminosity function and the statistics of beaming, 
and find that our simulated data produce results which are in¬ 
distinguishable from those produced by the real data on which 
they are based. 


Acknowledgements. It is a pleasure to thank Katherine Blundell for 
reading the manuscript prior to submission and providing insightful 
comments. 


3.5. Discussion 

Table presents the results of correlation and regression anal¬ 
ysis on the simulated datasets, and clearly shows that the sim¬ 
ulated datasets display exactly the same correlations as the real 
data. Whether it is core or extended radio luminosity that is 
correlated with core dominance is therefore not due to some 
physical difference between CDQs and LDQs, but merely an 
effect caused by the sample selection criteria, as discussed in 
We note that our simulated data underproduces the num¬ 
ber of highly core-dominated (i?^10®) quasars, compared to 
observations. This is presumably a result of our not having in- 
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Fig. 1. Plots of extended and core radio luminosity (integrated fr om 10 MHz-250 GHz) vs the ratio of these quantities, for a 
simulated dataset with the same selection criteria as that of Punsly (1995). 
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Fig. 2. Plots of extended and core radio luminosity (at 1.64 GHz in the rest fra me of the source), vs the ratio of these quantities 
for a simulated dataset with the same selection criteria as that of Murphy et al. (1992). 
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Fig. 3. Plots of extended and core radio luminosity (at rest-frame 5 GHz) vs the ratio of these quantities for a simulated dataset 
with the same selection criteria as that of Hough & Readhead (1989). 


Hough, D.H., Readhead, A.C.S. 1989, AJ, 98, 1208 
Isobe, T., Feigelson, E.D., Nelson, P.I., 1986, ApJ, 306, 490 
Laing, R.A., Riley, J.M., Longair, M.S. 1983, MNRAS, 204, 151 
Morganti, R., Oosterloo, T.A., Reynolds, J.E., Tadhunter, C.N., Mi- 
genes, V., 1997, MNRAS, 284, 541 
Murphy, D.W., Browne, I.W.A., Perley, R.A. 1992, MNRAS, 264, 298 
Orr, M.J.L., Browne, I.W.A. 1982, MNRAS, 200, 1067 
Punsly, B. 1995, AJ, 109, 1555 


Qin, Y.P., Xie, G.Z., Bai, J.M., Fan, J.H., Zheng, X.T., Wu, S.M. 1998 
, A&A, 333, 790 (Q98) 

Simpson, C., 1996, VA, 40, 57 
Simpson, C., et ah, 1998, in preparation 
Urry, C.M., Padovani, R, 1995, PASP, 107, 803 











